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(Dated: February 7, 2008) 

We present measurements of the branching fractions B(B° — > D*~ D* + ) and B(Df — > 4>tt + ), based 
on 123 x 10 6 T (4S) -> BB decays collected by the BABAR detector at the PEP-II asymmetric-energy 
e e~ B factory. A partial reconstruction technique is used to measure B(B° -> D*~D* S + ) and the 
decay chain is fully reconstructed to measure the branching fraction product B(B° —> D*~ D* + ) x 
B(D^ — + 4>tv + ). Comparing these two measurements provides a model-independent determination 
of the D+ -> (j}-K+ branching fraction. We obtain B{B° -> D*~D*+) = (1.88 ± 0.09 ± 0.17)% and 
B(D+ -> 07T+) = (4.81 ± 0.52 ± 0.38)%, where the first uncertainties are statistical and the second 
systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Published measurements of B(B a -> D*~ D*+) 0, 
are limited by the uncertainties on the D+ partial de- 
cay rates. A substantial improvement can therefore be 
obtained using a partial reconstruction technique where 
the Df is not explicitly reconstructed. The measurement 
of B(B° — > D*~D* + ) provides a test of the details of 
the factorization assumption Q in the relatively high q 2 
regime 4]. Partial reconstruction in addition allows an 
unbiased measurement of the Df — > (f>n + branching frac- 
tion, which has important implications for a wide range 
of D s and B physics, as most of the D s decay branching 
fractions are normalized to it PJ. As an example, an im- 
proved measurement of B(Df — > </>7r + ) would reduce the 
experimental uncertainty on the constraint on the Uni- 
tary Triangle parameter 7 from the measurement of the 
CP violating asymmetry in B° D^ttF decays 0. 

We used (123 ± 1) X 10 6 BB decays collected at 
the PEP-II asymmetric-energy e + e~ B factory with 
the BABAR detector, which is described in detail 
elsewhere . We provide here a brief description of the 
detector components relevant for this analysis. Charged- 
particle trajectories are measured by a silicon vertex 
tracker (SVT) and a drift chamber (DCH) immersed in 
a 1.5 T solenoidal magnetic field. The five-layer SVT 
enables tracks with low transverse momentum to be re- 
constructed. The energy and direction of photons and 
electrons are measured by a CsI(Tl)-crystal electromag- 
netic calorimeter (EMC). Charged-particle identification 
is obtained from the measurement of energy loss in the 
tracking system, and from the measurement of the num- 
ber and the angle of Cherenkov photons in a ring-imaging 
Cherenkov detector (DIRC). 

To study efficiencies and backgrounds and to validate 
the analysis we use several event samples produced with 
a Monte Carlo (MC) simulation of the BABAR detector 
based on GEANT4 [7| and reconstructed through the 
same chain as the data. 

The B° D*-D*+ -> {D+-f)(D°TT-) decay is 
reconstructed using two different methods. The first 
method combines the fully reconstructed D*~ decay with 
the photon from the D* + — > D+7 decay, without explicit 
-D+ reconstruction. Denoting the measured yield by Md 3 , 



we can write: 



B{B° -> D*-D*+) =Bi = 



(i) 



Here JC = 2N BS f oB(D*+ D+ 7 )B(D*~ -► D^r"), 
N B g is the number of B-meson pairs, /oo = 0.499 ± 
0.012 is the fraction of Y(4S) -> B°B° decays, B t is 
the branching fraction of D° decay mode i, Si is the effi- 
ciency for partially reconstructing the B° with a photon, 
a low momentum ("soft") pion and a D° reconstructed 
in mode i. 

The second method, based on full reconstruction of 
the B° -> D*-D* s + decay via D+ -> 0tt+ (0 K+K~), 
measures the branching fraction product B2 = B(B Q — » 
D*-D*+) x B(Df -> ^+): 



B 2 = 



Ml 



KB{4>^K+K-)YMBiY 



(2) 



where Nd 3 ->4>-k is the number of reconstructed decays 
and £^ is the efficiency for fully reconstructing the B°, 
including reconstruction of <j> — > K + K~ . The D+ — » 
4>ir + branching fraction is measured from the B2/B1 ratio: 



B2 
Bi 



(3) 

where the factor K, drops out. Although the efficiencies 
£j and are in general different, they include common 
factors and many systematic uncertainties cancel in the 
ratio. 

To extract the signal in partially reconstructed events, 
we compute the "missing mass" recoiling against the 
D*~j system, assuming that a B° — > D*~jX decay took 
place: 

m-miss = J {E B - E D * - E^) 2 - (p B -p D . -p 7 ) 2 , (4) 

where all quantities are defined in the T(4S) center-of- 
mass (CM) frame. While the photon and D*~ energies 
(Ej, Ed*) and their three-momenta (p 7 , Pd*) are mea- 
sured, kinematical constraints are needed to determine 
the B four- momentum (Eb, Pb)- In order to do that we 
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equate the B-meson energy with Sbeami the beam en- 
ergy in the CM frame, and calculate the cosine of the 
opening angle i^bd* between the B and the D*~ mo- 
mentum vectors from 4-momentum conservation in the 
B° — > D*~ D* s + decay. This leaves the azimuthal angle 
of the B meson around the D* ~ direction as the only 
undetermined parameter in the kinematics of the decay. 
MC studies show that an arbitrary choice of this angle 
(we fix cos </>bd* = 0) introduces a negligible spread (of 
the order of 1.5 MeV/c 2 ) in the m m ; ss distribution. The 
m m j ss distribution of signal events peaks at the nominal 
D+ mass [H with a width of about 15MeV/c 2 . 

We suppress unphysical D*~j combinations by requir- 
ing | cos#BD. | < 1.2 and events from e + e~ ^uu, dd, ss, 
cc production by requiring the ratio of the second to the 
zcroth Fox- Wolfram moments 10] to be less than 0.3. 



D* candidates are reconstructed in the D 



mode 



using D° decays to K + tt~ , K + ir~ir + ir~ , K + tt~tt° and 
K®ir + TT~ , listed here in order of decreasing purity. The 
X 2 probabilities of both the D° and D* vertex fits are 
required to be greater than 1%. The D*~ momentum in 
the T{AS) frame must satisfy 1.4 < p D , < 1.9 GeV/c. We 
require the reconstructed mass of the D° to be within 3 
standard deviations (<x m ) °f the measured peak value, 
andQ^* = mr>* — m D o — m OT to satisfy Q lo < Q D * < Q hi , 
where the choice of limits Q Xo = 4.10 - 5.20MeV/c 2 
and Qhi = 6.80 — 7.90MeV/c 2 around the nominal 
value Q™ G = 5.851 McV/c 2 depends on the D° decay 
mode. Kaon identification is required in K + ir~ir° and 
K + ir~ir + Tr~ modes. The from the K®ir + ir~ mode 
must have an invariant mass within 15MeV/c 2 of the 
nominal K® mass and a flight length greater than 3 mm. 

If more than one D*~ candidate is found, we first re- 
tain those that have the D° reconstructed in the decay 
mode with the highest expected purity. If ambiguities 
persist at this stage, we choose the best candidate based 
on the track quality of the soft pion and finally on the 



minimum value of x 2 = [(Qd* -Qd* )/ a Q D *] +[( m D° 



,PDG 
l D o 



where 



is the measured resolution 



on Qd'- 

Photon candidates are chosen from clusters of energy 
deposited in the EMC that are not associated with any 
charged track. The energy spectrum of photons from the 
D* s + -> D+7 decay is rather soft {E 1 < 0.4 GeV) and this 
makes controlling the background due to random photon 
associations one of the main challenges in the analysis. 
We require E 1 > 142 MeV and use the energy profile 
of the cluster to refine the photon selection, requiring 
a minimum cluster lateral moment |lll | of 0.016, and a 
minimum Zernike moment A20 ^2 °f 0-82. We also re- 
ject photon candidates that form in combination with 
any other photon in the event a 7r° whose invariant mass 
is between 115 and 155 MeV/c 2 and whose momentum in 
the CM frame is greater than 200 MeV/ c. This selection 
retains more than one photon candidate in about 10 % of 
the events. In these occurrences we choose the one that 



maximizes the value of a likelihood ratio based on the 
energy and the shape of the reconstructed cluster. 

The cuts are chosen to maximize the expected statis- 
tical significance of the selected signal using MC. The 
combinatorial background is dominated by B°B° events. 
None of the background components peak at the D+ 
mass in the m m i BS distribution. The reconstruction and 
selection efficiency, evaluated on simulated events, is 
(eB> = E 4 (oBi) = (5.15 ± 0.03) x 10~ 3 . 

We extract the signal yield using an unbinned 
maximum-likelihood fit to the m m i ss distribution. The 
signal peak is well described by a Gaussian proba- 
bility density function (p.d.f.). We parameterize the 
combinatorial background with the threshold function 



B(m miss ) = B (l 



j-(m max -m misB 



^ ) (^miss/^max) 



Fig. n shows the result of the fit to the missing-mass 
distribution. The width of the Gaussian signal distribu- 
tion is taken from MC simulation. The signal yield is 




2.05 

2 

Missing mass (GeV/c ) 

FIG. 1: Fit (solid line) to the measured missing-mass distribu- 
tion. The background component is shown as the dashed line. 

./Vd, = 7488 ± 342 events, corresponding to a branching 
fraction B(B° -> D*~D*+) = (1.88 ± 0.09) %, where the 
quoted error is purely statistical. 

We now describe the full reconstruction of the B° — > 
D*-D*+ -> (D+~/)(D°ir-) chain, with D° decaying 
into the four modes considered, and Df — > 
K + K~n + . Two kinematical variables are used: AE = 
Eb — -Eboam and the energy-substituted mass tties = 



V ^bcam — Ps • T ne t w0 variables have very little cor- 
relation; for signal events AE peaks around zero and 
mES at the B-meson mass. After applying selection cuts 
(described below) on the D* + and D*~ candidates, we 
retain the combination with the smallest value of \AE\. 
The number of fully reconstructed B° candidates is then 
obtained from a fit to the tties spectrum. 

The selection of D*~ candidates and most of the re- 
quirements on photon candidates are identical to those 
adopted in the partial reconstruction analysis. Due to the 
additional kinematical constraints on fully reconstructed 
B decays, the combinatorial background level is much 
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smaller; we can therefore relax the requirement on E~, 
thus improving the statistical significance of our sam- 
ple. We reconstruct <j> candidates from pairs of oppositely 
charged tracks, with at least one track satisfying kaon 
selection criteria; candidates are formed by combi- 
nation with an additional track, with charge opposite 
to the soft pion from the D*~ decay. A mass within 
±50 MeV/c 2 of the nominal D+ mass Q is required. Fi- 
nally, D*~ and D* + mass constraints are imposed in or- 
der to improve the toes and AE resolution of the B° 
candidate. We require the tod* — to,£> s mass difference 
to be between 125 and 160 MeV/c 2 , the reconstructed 
4> mass to be between 1.008 and 1.035 GeV/c 2 , EL to be 
greater than 90MeV, and |AE| to be less than 50MeV. 




5.2 5.22 5.24 5.26 5.28 5.3 

m ES (GeV/c 2 ) 



FIG. 2: Fit (solid line) to the measured m^s distribution. 
The background component is shown as the dashed line. 

We perform an unbinned maximum-likelihood ht to the 
Toes distribution with the sum of a Crystal Ball |l3| func- 
tion, and a threshold ARGUS |lj] function; the latter 
accounts for the combinatorial background. From the ht 
to the data sample, shown in Fig.0 we obtain (247 ± 19) 
events in the signal region defined as toes > 5.27 GeV/c 2 . 

MC studies indicate a peaking contribution due to real 
B° — > D*~D* + events, where either the D° does not de- 
cay into the reconstructed modes, or the Df does not 
decay into 4>tt + . We subtract the peaking background ap- 
plying a correction factor to take into account that the 
values of the B° -> D*~D*+ and D+ -> 0tt+ branch- 
ing fractions that we have measured are different from 
those used in the simulation, with an iterative procedure. 
The resulting number of peaking background events ex- 
pected in the data sample is 35 ± 6 events; this uncer- 
tainty is taken into account in the systematic error. Af- 
ter subtraction of the peaking background events, the 
final signal event yield is Afu s ^. ( j )n = (212 ± 19). Tak- 
ing into account the reconstruction and selection effi- 
ciency (e'B) = (4Bi) = (6-16 ± 0.24) x 10~ 3 , eval- 
uated on simulated events, we determine £>2 = B(B° — > 
D*-D*+) x B(D+ ^0tt+) = (8.81 ±0.86) x 10- 4 , where 
the error is statistical only. 



TABLE I: Summary of systematic uncertainties. 



Source 


Bi [%\ 


B 2 [%\ 


Ba/Bi [%] 


p.d.f. modeling 


4.8 




4.8 


Comb, background 




2.9 


2.9 


MC statistics 


0.6 


3.2 


3.3 


Peaking background 




2.8 


2.8 


B counting 


1.1 


1.1 




/oo 


2.4 


2.4 




Soft pion efficiency 


2.2 


2.2 




D*~ Tracking efficiency 


2.4 


2.4 




D*~ Vertexing efficiency 


2.0 


2.0 




-D+ Tracking efficiency 




2.6 


2.6 


D£ Vertexing efficiency 




2.0 


2.0 


Photon efficiency 


1.8 


1.8 




7r u eff. (D" -» K + n ty u ) 


1.2 


1.2 




7r veto 


4.7 


4.7 




Po rj--\ pip iH put i "H pai"i r»n 

J. CIL L 1L1C 1VJ.C11L111L/CILHJ11 


0.4 


0.4 




Polarization uncertainty 


0.8 


0.8 




D° branch, fract. [1] 


3.2 


3.2 




B(D*~ -> D°tt-) |4J 


0.7 


0.7 




B{D*+ - D+7) [15] 


0.8 


0.8 




B{(t>-> K + K~) [VJ 




1.2 


1.2 


Total systematic error 


9.1 


10.7 


7.9 



The main sources of systematic uncertainties on the 
B° — > D*~ D* + branching fraction measurement are 
listed in the second column (£>i) of Tabled We compared 
the resolution of the Gaussian p.d.f. in data and MC by 
fitting the missing mass distribution in the very clean 
sample of fully reconstructed B° — > D*~ D* s + events. We 
disentangle in this way the effect of the experimental res- 
olution on the width of the signal peak from the corre- 
lations in the fit between the width and the background 
parameters. We obtain adata/cMC = (1.01 ± 0.05). We 
repeated the TO m i ss fits changing the Gaussian width by 
this uncertainty, and varying the background parameters 
by their errors. We assign the maximum deviation as sys- 
tematic uncertainty, labelled in Table [I] as "p.d.f. mod- 
eling". The MC statistics uncertainty is the statistical 
error on the efficiency determination. The systematic 
uncertainties due to tracking, vertexing, photon and 7T° 
reconstruction efficiencies, and particle identification are 
evaluated using independent control samples. The ef- 
fect of the 7T° veto is evaluated from fully reconstructed 
events. The uncertainty due to the dependence of the ef- 
ficiency on the polarization of the B° — > D*~ D* + decay 
is assessed from MC samples generated with complete 
longitudinal and transverse polarization. In the full re- 
construction analysis the error on peaking background 
is due to the MC statistics and to the uncertainty on 
the relevant D° and Df branching fractions; the un- 
certainty on the combinatorial background is estimated 
using the AE sideband (|AE| > 200 MeV) as an alterna- 
tive way of computing the number of background events 
under the toes peak. Several systematic uncertainties in 
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the full reconstruction are in common with the partial 
reconstruction analysis, and therefore cancel in the ra- 
tio of Eq. |21 All remaining sources are listed in the last 
column of Table [I] 

We repeated both the partial and the full reconstruc- 
tion analyses on generic MC samples consisting of B°B°, 
B + B~ , and low-mass qq events, finding no bias. The 
result is also stable over different data-taking periods. 

In summary, we have measured the B° — > D*~D* + 
branching fraction 

B{B° -> D*-D*+) = (1.88 ± 0.09 ± 0.17)%, (5) 

where the first uncertainty is statistical and the second 
is systematic. This result is independent of the par- 
tial decay rates of the Df mesons. It is consistent 
with a previous BABAR measurement Q and with the 
world average, and reduces the total uncertainty by a 
factor of about three. The measurement is in agree- 
ment with the predictions of the factorization model 
B(B a -> D*-D*+) thcol = (2.4 ± 0.7) % 0. 

We have measured the branching fraction of Df — > 
<ftir + decay: 

B(Df -> <jm + ) = (4.81 ± 0.52 ± 0.38)%. (6) 

This result represents an improvement by about a factor 
of two over previous measurements |]J, [ljj . 
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